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COMPARISON OF REMOTELY SENSED CONTINENTAL-SHELF WAVE 

SPECTRA WITH SPECTRA COMPUTED BY USING A 

WAVE REFRACTION COMPUTER MODEL 

Lament R. Poole 
Langley Research Center 


SUMMARY 

An initial attempt to verify the Langley Research Center and Virginia 
Institute of Marine Science mid-Atlantic continental-shelf wave refraction 
model was made by using the model to simulate refraction occurring during a 
continental-shelf remote sensing experiment conducted on August 17, 1973- 
Simulated wave spectra compared favorably, in an overall qualitative sense, 
with experimental spectra. It was observed, however, that most of the wave 
energy resided at frequencies higher than those for which refraction and 
shoaling effects were predicted. In addition, variations among the experi- 
mental spectra were so small that they were considered statistically insig- 
nificant and, thus, not attributable to refraction and shoaling effects. In 
order to verify the refraction model, simulation must be performed in conjunc- 
tion with a set of experimental spectra in which significant variation among 
the individual spectra is exhibited and in which a considerable portion of 
the total energy resides at frequencies for which refraction and shoaling 
effects are expected . 


INTRODUCTION 

Synoptic monitoring of ocean surface conditions from the NASA SEASAT-A 
satellite (ref. 1) promises to have a significant impact on future coastal- 
zone planning and management activities. Analytical models under development 
at the present time will use the broad-scale data from Seasat to provide 
short-term predictions of ocean surface conditions on a smaller spatial 
scale. These models could be of aid to the coastal-zone manager in areas 
such as disaster warning, pollution dispersal, or site selection for offshore 
and coastal facilities. 

The wave refraction model is one analytical procedure which could play 
an integral role in an effective coastal-zone management program. Given ref- 
erence wave measurements at points along the outer edge of the continental 
shelf, the refraction model would predict the behavior of the waves as they 
cross the shelf and impinge on the shoreline. Of course, before a wave 
refraction model can be used in an operational management or planning pro- 
gram, its computations must be verified by comparison with actual continental- 
shelf wave measurements for a wide range of wave conditions. 


ii ■■ it i i i ill ii i ii ii ii i a i ii i 1 1 1 ■ i i I ll III l 


I 



The present paper describes a procedure used in an initial attempt to 
verify the Langley Research Center (LaRC) and Virginia Institute of Marine 
Science (VIMS) wave refraction model described in reference 2. Experimental 
wave measurements were obtained by using an airborne laser profilometer 
(ref. 3) which was flown over an approximately 60-nautical-mile stretch of 
the continental shelf off Assateague Island, Maryland, on August 17, 1973. 
Energy density spectra derived from the experimental measurements are com- 
pared with simulated spectra computed by using the LaRC-VIMS refraction 
model according to the procedure of Chao (ref. 4). Discussion of agreement 
between computed and experimental spectra leads to recommendations for future 
attempts at model verification. 


SYMBOLS 

b ray separation distance, meters 

d water depth, meters 

E energy density between ray pair j and j + 1 , joules/meter 2 

F(f) normalized energy amplification function 

f wave frequency, hertz 

H.j /3 significant wave height, average height of highest one-third of 

waves, meters 

j wave ray index 

K r refraction coefficient, defined by equation (2) 

K g shoaling coefficient, defined by equation (3) 

k wave number, meters - ^ 

S(f) spectral density, meters^/hertz 

T wave period , seconds 

x,y distances along axes of wave refraction model, nautical miles 

a angle from which waves are propagating, as defined in figure 4, 

degrees 

Subscript : 

o initial (deep-water) value 

A bar over a symbol indicates spatial average over a flight-track 
segment . 
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WAVE MEASUREMENT EXPERIMENT 


The nominal plan chosen for airborne continental-shelf wave measurement 
experiments consists of two phases. During the first phase the aircraft 
flies offshore for a distance of approximately 60 nautical miles in a direc- 
tion opposite to that of the local (onshore) wind and at an altitude of 
1.6 kilometers. Along this flight track, reference photographs are taken of 
the ocean surface with an aerial mapping camera. The second phase consists 
of a return (onshore) flight in the direction of the local wind and at an 
altitude of 160 meters. During this phase a continuous profile of the ocean 
surface along the flight track is obtained by using a laser profilometer in 
the mode described in reference 3* 

The present wave measurement experiment was conducted on August 17, 

1973 > aboard a NASA 0-5** aircraft stationed at Wallops Flight Center, Vir- 
ginia. Weather conditions in effect during and for several hours in advance 
of the flight included scattered clouds at an altitude of 650 meters and a 
visibility of 6 to 7 nautical miles. Winds were from the northeast, shifting 
to the east-northeast at speeds of 5 to 8 meters/second. On the basis of 
these wind data, the nominal bearing for the offshore flight track was 
selected as 070° from true north. A navigational buoy located off Assateague 
Island, Maryland, near longitude 75° W and latitude 38° N was chosen arbi- 
trarily as a geographical reference point over which to fly and initiate the 
experiment . 

The offshore phase of the experiment was begun at 17:**6 GMT on August 17. 
Onboard navigation was provided by a L0RAN navigation system with which read- 
ings were made at 10-nautical-mile intervals to allow for course corrections 
as necessary. Onboard L0RAN readings taken at points from 0 to 60 nautical 
miles offshore from the reference buoy were mapped onto a Mercator projection 
by using a U.S. Naval Oceanographic Office L0RAN navigational chart. Shown 
in figure 1 is the offshore flight track faired through the actual data 
points. Some deviation from the nominal 070° bearing is evident in figure 1. 
For the purposes of this paper, the possibility of errors in the L0RAN read- 
ings is disregarded and the faired flight track is considered exact. Refer- 
ence photographs were taken during this offshore flight with the aerial map- 
ping camera equipped with a yellow filter and 9-inch format film. Three of 
the reference photographs are presented in figure 2, with parts (a), (b) , 
and (c) first (shoreward), middle, and last portions of the flight track, 
respectively. In each figure, the aircraft was flying in a direction from 
left to right (offshore), and the majority of waves are traveling in nearly 
the same direction . 

The offshore phase of the experiment was concluded at 18:06 GMT at a 
point approximately 60 nautical miles from the reference buoy. The aircraft 
was then brought down to the 160-meter altitude for the return flight, which 
was begun at 18:15 GMT at a nominal bearing of 250° from true north and a 
nominal speed of 80 meters/ second . Onboard LORAN readings taken at 10- 
nautical-mile intervals on the return flight were again mapped onto a Merca- 
tor projection. The actual onshore flight track faired through the projected 
LORAN data points is shown in figure 3. Considerable deviation from the nom- 
inal bearing is evident, but, again, the faired flight track is considered 
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exact. During this phase of the experiment, a continuous profile of the 
ocean surface (with a maximum resolution of 10 cm) was obtained by using the 
laser profilometer . The laser data were recorded in analog form on magnetic 
tape between the times of 18:15:30 GMT and 18:35:30 GMT. 


LASER DATA ANALYSIS 

In order to provide quasi-homogeneous data sets for the subsequent spec- 
tral analysis, the magnetic tape record of laser profilometer data was parti- 
tioned into 10 adjacent but independent sections of 2-minute length each. 

Each section was digitized using an interval of 0.04 second and then treated 
by a high-pass digital filter (ref. 5) in order to remove extraneous energy 
at very low true frequencies (less than 0.06 Hz) attributed to vertical air- 
craft motions. On the basis of evidence from a similar wave measurement 
experiment (ref. 6), aircraft pitch and roll motions were assumed to have lit- 
tle influence on the data record obtained. Raw apparent spectra (spectra 
observed from the moving aircraft reference frame) for each filtered section 
were obtained by computing the Fourier cosine transform of the autocorrela- 
tion function (to 150 lags) of the filtered section in question, according 
to the method of Blackman and Tukey (ref. 7). Such a treatment of data 
obtained over a finite spatial increment produces, in essence, a spectrum 
which is a spatial average of local spectra over the increment in question 
(ref. 8). The raw apparent spectra were then treated by the Hamming weighted 
averaging process (ref. 9, pp. 447-448) to obtain smoothed apparent spectra 
with 80 percent confidence that repeated measurements of the same wave field 
would result in spectral density values lying between 0.71 and 1.31 times the 
values computed from the initial measurements. 

For the purposes of the present paper, it was assumed that directional 
spreading in the wave field and deviations of the aircraft flight track from 
the mean direction of wave propagation could be neglected. Therefore, the 
smoothed apparent spectra were transformed into true spectra in a stationary 
reference frame by using the unidirectional transformation technique pre- 
sented in reference 10. In order to obtain spatially averaged depth values 
along the various flight-track segments, which were required for the trans- 
formation, the onshore flight track was mapped into the rectangular coordi- 
nate system of the LaRC-VIMS wave refraction model (fig. 4). The flight 
track was divided by tick marks into 2-minute segments, which are numbered 
from 1 to 10 (advancing shoreward), corresponding to the previously computed 
apparent spectra. Segment 1 was found to lie outside the seaward boundary 
of the refraction model coordinate system and was selected to provide the 
reference spectrum to be used in subsequent calculations. The average depth 
d along segment 1 was estimated from an auxiliary bathymetry chart to be 
100 meters, and the resulting transformed spectrum, designated as the refer- 
ence deep-water spectrum, is shown in figure 5- 

A portion of segment 2 was also found to lie outside the seaward bound- 
ary of the model coordinate system. Since there would be no capability for 
simulating in full the propagation and refraction of the reference spectrum 
through segment 2, the apparent spectrum for this segment was discarded. 
Average depth values along segments 3 to 10 were computed by averaging depth 
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values at points along appropriate rays in the wave refraction diagrams which 
were constructed for use in simulating refraction of the reference deep-water 
spectrum (to be discussed in the following section). The resulting trans- 
formed spectra for segments 3 to 10 are shown in figure 6. 


SIMULATION OF SPECTRUM REFRACTION 

Because of the weak local winds in effect during the onshore flight 
experiment, it was assumed that any changes observed in the spectra obtained 
over the flight track were due to refraction and shoaling effects rather 
than to local wave generation. Simulation of these refraction and shoaling 
effects was performed by using the LaRC-VIMS monochromatic wave refraction 
model according to a procedure such as the one outlined by Pierson, Neumann, 
and James (ref. 11) and implemented by Chao (ref. 4). It was assumed that 
the deep-water spectrum could be represented, with wave-wave interaction 
neglected , by a linear superposition of monochromatic waves traveling ini- 
tially in a single direction. (This assumption correlates with the neglect 
of directional spreading in analysis of the laser data.) The first step in 
the simulation procedure was then to construct a refraction diagram for each 
monochromatic wave used to approximate the spectrum . 

The LaRC-VIMS refraction program was altered for convenience in construc- 
tion of the required refraction diagrams. Computations were limited to a 
family of ' wave rays (orthogonals to wave crests), each separated by 1 nauti- 
cal mile in deep water, in the vicinity of the onshore flight track. Instruc- 
tions were added to the program to initiate ray computations at deep-water 
points lying along a crest line orthogonal to the initial ray direction in 
order to reference all computations to an intial time. Succeeding computa- 
tions were then performed in fixed time increments, and tick marks were added 
to rays in the constructed refraction diagram at spatial intervals correspond- 
ing to 5 time steps. The initial angle ct Q selected for all waves in the 
approximating series was the mean acute angle between the onshore flight track 
and the X-axis of the model coordinate system, which was computed to be 67°. 
Refraction diagrams were then computed for waves propagating from this direc- 
tion and ranging in period from T = 6 seconds to T = 14 seconds. This 
period range corresponds to a frequency range from f = 0.1667 hertz to 
f = 0.0714 hertz. Little energy was observed in the transformed experimental 
spectra (figs. 5 and 6) at frequencies lower than 0.0714 hertz, and waves at 
frequencies higher than 0.1667 hertz were found to undergo little or no 
refraction. Sample refraction diagrams are presented in figure 7 for periods 
T = 6, 8, 10, 12, and 14 seconds. The circular symbols in each figure denote 
LORAN aircraft position data for the onshore flight track . 


Energy Amplification Functions 

Variations in wave energy due to refraction and shoaling effects can be 
estimated according to linear theory by considering the wave energy between 
two vertical planes which are orthogonal to the wave crests and intersect 
with the surface to produce wave rays. If reflection, percolation, and bot- 
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tom friction are neglected and it is assumed that no energy is transmitted 
across wave rays, then for a wave of frequency f n , the energy density E 
at some point shoreward from a reference deep-water location and bounded by 
ray pair j and j + 1 can be written as 

E ( f n) = W K r 2 ( f n> K S 2 ( f n> (1) 

where E Q is the initial (deep-water) energy density between ray pair j 
and j + 1; K r is the local refraction coefficient, defined by 



where b Q is the initial (deep-water) distance between rays j and j + 1 
(fixed as 1 n. mi.) and b is the distance between rays j and j + 1 at 
the local point of interest; and K g is the local shoaling coefficient, 
defined by (ref. 2, p. 8) 

/ ? \ 1/2 

K = | 2 cosh g kd \ (o 

s \2kd + sinh 2kd/ 


where d is the depth at the point of interest . By computing the energy 
density ratios E/E., at the point of interest for waves of varying fre- 
quency, a normalized energy amplification function F(f) can be constructed 
as a function of wave frequency. Then, with a reference deep-water spectral 
density S Q (f) given, the local unidirectional refracted spectral density 
S(f) can be computed as the product of the reference spectral density and 
the normalized energy amplification function, or, 


S(f) = S Q (f) F(f ) 


(4) 


However , the wave spectra computed from the experimental laser profiles 
represent spatial averages of local spectra over the various flight-track 
segments, as discu ssed previously. Accordingly, spatially averaged amplifi- 
cation functions F(f) were computed for use in the simulation procedure by 
first superimposing the aircraft flight track onto each of the constructed 
refraction diagrams as shown, for example, in figure 8 for a wa ve p eriod 
T = 10 seconds. A spatially averaged normalized energy ratio E/E 0 for seg- 
ment 8, for example, was computed by averaging local energy ratios computed 
by the refraction program at points between p 1 and p 2 * Performance of 
such computations for each monochromatic wave used to approximate the spec- 
trum (T ranging from 6 to 14 sec) resulted in construction of average energy 
amplification functions for each segment, in the frequency range from 0.0714 
hertz to 0.1667 hertz. For frequencies less than 0.0714 hertz, the respec- 
tive amplification functions were assumed to be equal to their values at 
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0.0714 hertz, while for frequencies greater than 0.1667 hertz, the amplifi- 
cation functions were assumed to be 1.0. Plots of the amplification func- 
tions for the eight flight-track segments are given in figure 9, in which the 
dashed portions of the plots indicate assumed values. 


Simulated Refracted Spectra 

Simulated refracted spectra for flight-track segments 3 to 10 were com- 
puted by multiplying, according to equation (4), the reference deep-water 
spec trum (fig. 5) by the spatially averaged energy amplification function 
F(f) for each of the segments. The resultant simulated spectra are compared 
with the experimentally determined spectra in figure 10. 


DISCUSSION 

In an overall qualitative manner the simulated spectra for all eight 
flight-track segments compare favorably with the respective experimental spec- 
tra. With respect to verification of the refraction model, however, the 
results are less favorable. It can be seen that nearly all the energy in the 
various spectra resides at frequencies greater than 0.12 hertz. Since the 
computed amplification functions were essentially 1.0 in this frequency range 
for all flight-track segments, the simulated spectra are identically equal to 
the reference deep-water spectrum in this range. In addition, if the spectra 
in figure 10 are compared with the reference deep-water spectrum (fig. 5), 
the spectra for each flight-track segment appear to lie, for the most part, 
within the 80-percent confidence limits of the reference spectrum. Thus, in 
an overall manner the changes in the various spectra relative to the refer- 
ence deep-water spectrum are statistically insignificant at the 80-percent 
confidence level . 

Since the various spectra exhibit a similar characteristic shape, dif- 
ferences among the spectra can be discussed from a more quantitative approach 
by considering the significant height under the assumption of 

a Gaussian distribution of wave heights, is equal to 4 times the square root 
of the area under the spectrum (ref. 12, pp. 22). The areas under simulated 
and experimental spectra for flight-track segments 3 to 10 were computed by 
trapezoidal integration, and the resultant significant heights were compared 
with the significant height for the reference deep-water spectrum. Relative 
changes of 0 to 2 percent in the simulated significant heights substantiate 
the observation that most of the energy in the reference spectrum resided at 
frequencies above those for which refraction effects were predicted. A maxi- 
mum relative change of 9 percent was found in the experimental significant 
heights. Such a change would correspond to a 19-percent maximum difference 
in spectral density levels, which falls within approximately +.30-percent 
band limits that must be exceeded to attain a statistically significant 
difference . 

The foregoing discussion indicates that in order to verify the LaRC-VIMS 
refraction model , experimental data must be obtained under conditions such 
that a considerable portion of the total wave energy resides at frequencies 
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lower than those experienced in the present analysis. In addition, there 
must be variations among the individual wave spectra to a degree that the 
variations can be attributed with statistical significance to physical 
causes, namely shoaling and refraction, rather than to random fluctuations. 

To achieve statistical significance while retaining reasonable frequency res- 
olution, future analyses might require much longer data record lengths or 
digitization using a much smaller interval . 


CONCLUDING REMARKS 

The present paper has described the procedure used in an initial attempt 
to verify the Langley Research Center and Virginia Institute of Marine Sci- 
ence mid-Atlantic continental-shelf wave refraction model. Experimental 
spectra were obtained along a 60-nautical-mile track in the region off 
Assateague Island, Maryland, on August 17, 1973, by using a laser profilom- 
eter mounted aboard a NASA C-5^ aircraft. The refraction model was then used 
in a linear superposition mode to simulate the spectrum refraction occurring 
during the flight experiment . 

The simulated spectra compared favorably with the experimental spectra 
in an overall qualitative sense. It was observed, however, that most of the 
wave energy resided at frequencies higher than those for which refraction 
effects were predicted. In addition, variations among the experimental spec- 
tra were statistically insignificant and thus could not be attributed to 
physical causes, namely shoaling and refraction. In order to verify the 
refraction model, simulation must be performed in conjunction with a set of 
experimental spectra in which significant variation among the individual 
spectra is exhibited and in which a considerable portion of the total wave 
energy occurs at lower frequencies than experienced in the present analysis. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
November 22, 1976 
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Figure 1.- Offshore flight track for wave measurement experiment. 
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Figure 2.- Continued 
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Figure 2.- Concluded. 
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Figure 3-- Onshore flight track for wave measurement experiment. 
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Figure 5.- Reference deep-water wave spectrum, transformed 
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Figure 6.- Wave spectra for flight- track segments, transformed for d 

of respective segments . 
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Figure 6.- Continued. 
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(c) Segment 5; d = 51 meters. 


Figure 6.- Continued. 
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(c) T = 10 seconds. 
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Figure 7.- Continued 
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(e) T = 14 seconds. 
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Figure 7.- Concluded 
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Figure 9.- Spatially averaged energy amplification functions for flight- 
track segments. Dashed curve indicates assumed values. 
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Figure 9.- Continued. 
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Figure 10.- Continued 
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(c) Segment 5. 


Figure 10.- Continued. 
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( f ) Segment 8 . 

Figure 10.- Continued. 
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(g) Segment 9. 
Figure 10.- Continued. 
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(h) Segment 10. 


Figure 10.- Concluded. 
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